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Abstract. We study an asymptotic behavior of the return probability for the critical 
random matrix ensemble in the regime of strong multifractality. The return probability 
is expected to show critical scaling in the limit of large time or large system size. Using 
the supersymmetric virial expansion we confirm the scaling law and find analytical 
expressions for the fractal dimension of the wave functions d,2 and the dynamical 
scaling exponent /z. By comparing them we verify the validity of the Chalker's ansatz 
for dynamical scaling. 
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1. Introduction 

It is well-known that the wave functions at the point of the Anderson metal-insulator 
transition are fractal [I]. Their amplitudes exhibit self-similar fluctuations at different 
spatial scales. The standard way to quantify such a complicated behavior is to consider 
the scaling of moments of the wave functions ip n {^) with the system size L: 

I q = Y,{\^n{v)\ 2q )^L- d ^- l \ (1) 
r 

where (...) stands for averaging over disorder realizations and over a small energy 
window. The fractal dimension d q , which is different from zero and from the 
dimensionality of the space d, is a fingerprint of the fractal wave functions. For the 
multifractal wave functions d q depends non-trivially on q, thus an infinite set of scaling 
exponents is required for the full description of the wave functions in this case. 

Additionally to the non-trivial scaling of the moments of the critical wave functions 
taken at a fixed energy, the correlations of wave functions at different energies show the 
critical behavior as well. The simplest correlation function involving two eigenstates 
corresponding to two different energies E m and E n can be defined as 

C(u>) = J2(\'Mr)\ 2 [K{r)\ 2 5{E m -E n - u>)). (2) 

r 

As any other correlator at criticality C (u) is expected to decay in a power-law fashion 

C(u>) oc {Eo/uY, A < u < E , (3) 

where A is the mean level spacing and Eq is a high-energy cutoff. What is more 
surprising is the fact that the dynamical exponent \i is related to the fractal dimension 
g?2 in a simple way 

H = 1 - d 2 /d. (4) 

This relation was suggested by Chalker and Daniel [21 [3] and confirmed by a great 
number of computer simulations (2j HJ E] thereafter. As Eq/uj 1 and /x > 0, Eq.(j3j) 
implies an enhancement of correlations in critical systems [I] which is possible only 
if there is a strong overlap of different wave functions. This rather counterintuitive 
picture becomes particularly striking in the regime of strong multifractality d 2 <C d, 
where eigenstates are very sparse [6] and almost localized. We are not aware of any 
analytical calculations supporting its validity in this case. 

The critical enhancement of the correlations plays an important role in the theory 
of interacting systems (cf. "multi-fractal" superconductivity [7j and the Kondo effect 
[8]) and, therefore, the theory of the critical correlations still attracts a considerable 
attention in spite of its long history. 

The aim of the present work is to demonstrate that the critical scaling holds true 
and confirm the validity of the dynamical scaling relation (J4j) with accuracy up to the 
second order of the perturbation theory. Some results of this work have been announced 
in a brief form in Ref . [12] . 
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The knowledge of the second order perturbative results is extremely important 
because it allows one to confirm the critical scaling. Besides, subleading terms in 
the scaling exponents can reveal some model-dependent features in contrast to the 
leading order perturbative result which is universal for a wide class of different critical 
models 0HU]. 

In this paper we consider a particular model relevant for the Anderson metal- 
insulator transition - the power law random banded matrix ensemble [11]. The matrix 
elements H mn of the Hamiltonian are given by the independent, Gaussian distributed 
random variables with the only constraint that matrix H is Hermitian. Their mean 
values are equal to zero and the variances are defined as 

2 1 J 1, n = m 

{lHmnl } = 2 \ \n - m\ » max{6, 1} (5) 

The long-range power-law decay of typical matrix elements \H mn \ 2 leads to the critical 
behavior at any value of the parameter b. This allows one to study the model 
perturbatively either for b ^> 1 or for b <C 1. The latter condition corresponds to 
the regime of strong multifractality investigated in the present work. In this regime the 
both scaling exponents /i and c?2 can be expanded into the power series in b and then 
compared term by term. 

The structure of the paper is as follows. In Section [2] we give an equivalent 
formulation of the Chalker's ansatz in terms of the return probability. In Section [3] the 
integral representation for the first order result for the return probability is derived. It 
is used in Section H] to calculate the scaling exponents in the first order in b. The second 
order result for the return probability is discussed in Section [5j The corresponding 
second order expressions for the dynamical exponent /i and fractal dimension d 2 are 
derived in Section [6] and Section [7] respectively. All calculations in Sections [3H7] are 
done for the unitary symmetry class. A brief announcement of analogous results for 
the orthogonal case is presented in Appendix B more detailed study will be published 
elsewhere. 



2. Return probability and scaling exponents 

It is convenient to reformulate the Chalker's ansatz (jlj) in terms of the return probability 

/oo 
due- iuJt C(co), (6) 
-oo 

where a matrix size N plays the role of the system size. Using the definition of C(u) 
([2]) it is easy to show that in the limit t — > oo the return probability tends to a finite 
limit, which is nothing else than the inverse participation ratio I2 

lim P N (t) = I 2 oc N~ d \ (7) 

Thus the knowledge of the return probability gives a way to calculate the fractal 
dimension g?2. On the other hand considering the limit iV — > 00 at a fixed large time t 
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one finds that 

lim P N (t) oc f- 1 , (8) 

iV/t-Kx> 

as it follows from Eq.((3]). So the dynamical scaling exponent \x can be also extracted 
from the behavior of the return probability. Thus we conclude that the Chalker's ansatz 
(]!]) is equivalent to the following statement 

d d 
/x-l=lim— — - lim hxP N (t) = lim lim In i>jv(t) = -d 2 . (9) 

t^oo 9 In t N/t^oo N^oo (9 In N t/N^oa 

In the regime of strong multifractality the return probability can be calculated 
perturbatively using the method of the virial expansion in the number of resonant states, 
each of them being localized at a certain site n. The virial expansion formalism was 
developed in Ref . [13] following the initial idea of Ref . [H] . The supersymmetric version 
of the virial expansion [TS] [TB] is formulated in terms of integrals over super-matrices. 
In particular, it allows us to represent P/v(t) as an infinite series of integrals over an 
increasing number of super-matrices associated with different sites: 

P N (t) = P (1) + P (2) + P (3) + . . . , (10) 

with PW = 1 and P^ = b 1 " 1 f^ l \bt). Functions f^' are governed by a hybridization of 
j localized states and can be calculated explicitly by means of integrals over j different 
supermatrices. The above expansion implies the corresponding expansion for lnP/v(i): 

\nP N (t) = pV+^- 1 -(P^+..., (11) 

where the first term is of the first order in b, two terms in the brackets are of the second 
order in b and so on. This representation allows one to find the corresponding power 
law expansions for the fractal dimensions g?2 and the dynamical exponent /i and hence 
to check the Chalker's ansatz in the form of Eq.Q order-by-order. 

We emphasize that the scaling exponents are finite constants and hence according to 
Eq-© the leading at N — > oo or t — > oo terms in lnP/v(t) must diverge logarithmically 
with only linear in IniV or \nt terms present. Higher order terms of the form ln m N 
or ln m t (m > 1) would generate divergent contributions to c?2 and fi indicating a 
violation of the power law scaling. As a matter of fact, both P^ and (p( 2 )) do 
contain higher order terms, such as In 2 N or In 2 1. We prove below that these divergent 
terms cancel out in the combination (p^ — | (p( 2 )) 2 j. This is a necessary condition 
for the existence of the critical scaling and of the Chalker's ansatz. 

3. Integral representation for P^ 

The return probability P(t) can be expressed in terms of the Green's functions as 

A 2 N 

P ^)^E«))). (12) 
p=i 
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where ((ab)) = (ab) — (a)(b) and the diagonal matrix elements of time dependent 
correlator Q pp (t) are related to its energy dependent counterparts Q pp (u) by the Fourier 
transform 

Gpp(t) = ^J due- tut Reg pp {uj). (13) 

For the latter quantity, defined by the product of the matrix elements of the retarded and 
the advance Green's functions Q pp {u) = G pp (E + u/2)G pp (E — u/2), the perturbation 
theory has been developed in Ref.[T5] . The leading order term of the perturbation 
theory, corresponding to the diagonal part of a random matrix, is KeQ pp (u) = 
(2-7T 2 / A)8(uj). Substituting this expression into Eq. f|T2|) yields = 1 reproducing the 
correct normalization of the return probability. The next order approximation taking 
into account an "interaction" between pairs of resonant states is given by Eq.(55) of 
Ref . [15] . The corresponding results for the return probability reads 

p(2) = 2 V^ {-2b pn t 2 ) k k 

Nt (k-lV. 2k -1 1 ' 

nj^p k=l v 

This expression was derived for an arbitrary variances of the off-diagonal matrix elements 
bpn = \{\H pn \ 2 )] for PLBRM model 6 pn = |(1 + \p - n\/b)~ 2 . In the large N limit, the 
double sum in Eq.( fl4|) may be replaced by the integral: 

N N fN ry pN 

VV/(b-n|)«2/ dy dxf(x)^2N dxf(x), (15) 



where the last equality is justified in Appendix A In the continuum limit, the 



counterpart of b pn is given by b 2 /Ax 2 (which is valid for \p — n\ ^> b), however this 
expression leads to the appearance of divergent integrals at x — > and hence should 
be regularized. To this end we replace b pn by 6 2 /4x 2 ^ 1_e - ) with e > and take the limit 
e — t- at the end of the calculations. Thus in the continuum limit we obtain 

, 2) r N , ^ / -b 2 t 2 \ k k i 



Now it is convenient to represent the last fraction as an integral 
Jp 1 dj3/3 2k ~ 2 , k > 1 and substitute this formula into Eq.f fT6|) 

-(3 2 b 2 t 2 \ k k 



i 



2fc-l 



2x 2(l-e) J (fc- 1)!' 



(17) 



Changing j3 by = f3bt/y/2 and using the fact that Y^k=i (~v) k (k-iy. = ~ v) e ~ 
we arrive at the following integral representation for the return probability: 



2 V 



Pm = [% l dxF2 (J^)i F 2 (y) = -^by 2 (l-y 
where r = bt/ y/2. 

| The right hand side of Eq.(55) should be multiplied by a/2, as in the present calculations we fix 
E = 0, while the averaging over E was performed in Ref.|15j 
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4. Scaling exponents: the first order perturbation theory results 

The above representation for P^(t) (118j) is a convenient starting point for calculation 
of d 2 and /i. The exponent [i can be extracted from the limit N — > oo of Eq. (fT8]) 



n 2 (r) = lim P( 2 ) = / f d yF 2 (AS) , y = (19) 

io Jo W J 



Differentiating Il2(r) with respect to lnr we obtain 
dU 2 



r—J, J= d yF 2 [ — \. (20) 



<91nr J \U 1 ~ e 

The last step in calculating /i is to take the limit e — > in the expression for Jy^j:- The r 
dependent factor then gives one and what we need to know is just the zeroth order, i.e. 
e independent term, in the e-expansion of J. The required expansion can be found with 
the help of the following general formula, which can be proved using the integration by 
parts: 

1 f°° df 

dpp 8 - 1 !^) = -/(0) - J rf/31n/3^ + 0(5). (21) 

To this end we change the integration variable y by t = l/y 2 ^ 1 "^ and apply the above 
formula: 

Vlirb r°, (1 -t)e~* irb r , 7 , x i ^, 9x , . 

J = ~ r/ tft^ {— = --= l + e(2 + ^ + ln2) + 0(e 2 ), (22) 

I 1 ~ e ) Jo t V2 L 2 J 

where 7 is the Euler's constant. We keep the first order in e term, as it is important in 

the second order perturbation theory. The zeroth order term gives the exponent fi: 

l-/i = ^ + 0(6 2 ). (23) 

Now we perform similar calculations for the fractal dimension d 2 - First we introduce 
in Eq. (ll8p the new integration variables x and j3 defined by the relations x = P'^Nx, 
(3 = N 1 ~ e [5 and then take the limit r — > 00 

1 

J dxF 2 ( —j . (24) 

Taking the derivative with respect to lniV and returning to the previous notation for 
the integration variables x — > x, (3 — > (3 we obtain: 

^""•r <25) 

Then we apply Eq. ([2Tj) with S = •jz - and find 



f/3 

d(3(3^ I dxF 2 



2e-l 



1 \ r°° ... „ 1 . 



dxF 2 ( —J + e jl rf/31n/3— F 2 (f3) + 0(e 2 ) ) . 
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The e-expansion of the first integral is given by Eq. f!22j) . while the second integral can 
be calculated explicitly using the definition of F 2 : 

J d (3\n(3-F 2 ^) = — (I + 7/2 + ln2). (27) 
Thus we arrive at the following result for 

<9n 2 



{21 



dlniV 

It is interesting to note that the first order in e term is absent in the above expression. 
The constant term yields the fractal dimension d 2 [32] 

d 2 = ^= + 0(b 2 ). (29) 

Comparing this result with the corresponding expression for 1 — fi Eq.( l2"B"j) we conclude 
that the Chalker's ansatz is valid in the first order perturbation theory. 

Leading contributions of order of b to the scaling exponents 1 — fi and d 2 in the 



orthogonal case are calculated in Appendix B 



5. Integral representations for 

The second order perturbation result for the return probability can be derived from the 
corresponding expression for the matrix elements of Green's functions given by Eq.(72) 
of Ref.[l5] : 

N N oo 

pm = ^Il E E (-8<w 2 ) fe (-«w 2 ) fa (-sw 2 )' 3 x 

p=l {m,n^p\ fci,2,3=0 

(*! + **)(*! + **- 1), (30) 



where 



Y{2[k x + k 2 + k 3 ] - I) 



*(k k k)-fl m - V2) 

"(«l.«2,«3j - 11 j/h.i 



W 2 fcj 
i=l 

(2hk 2 k 3 - hk 2 - hk 3 - k 2 k 3 ) T(h + k 2 + k 3 - 1). (31) 

First we multiply and divide the last expression by 2{k\ + k 2 + k 3 — 1) and then use the 
identity T{2z) = -±=2 2z - l / 2 T{z)T(z + 1/2) for z = k x + k 2 + k 3 - 1/2, which allows us 
to cancel T{ki + k 2 + k 3 ): 

00 

p(3) = F(k u k 2 ,k 3 )x 

fel,2,3=0 

(2*4 A;^ - fcife - fci&a - k 2 k 3 )(h + fc 2 )(fci + k 2 - 1) (32) 

with 

N N 

F(h, k 2 , k 3 ) = -J- ]T ^ (-86 pm t 2 ) fcl {-8b pn t 2 ) k2 (-%b mn t 2 ) ks x 



16t 2 iV 

p=l {m,n=£p\ 
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All terms in Eq. (l32j) are symmetric functions of k\, k 2 and k%, except for the term 
(ki + k 2 )(ki + k2 — 1). Symmetrizing it we can write P^ in the following form 

CO 

p(3) = J2 F (h,k 2 ,h) ^ A^aa^r^ 2 ^ 3 - (34) 

fcl, 2,3=0 Ol,2,3 

The coefficients A aiOC20C3 are invariant under permutations of the indices and hence all 
non-zero coefficients can be obtained from the following six: 

A)i2 = 2/3, A 013 = -2/3, A ni = 2, 

Am = -10/3, A 113 = 4/3, A 122 = 4/3. (35) 

The next step is to replace the summation over m, n, and p by the integration 
similarly to how it was done in the calculation of . To this end we notice, that 

f(\p — m\, \p — n\, \m — n\) « 

pN r-y r-y 

6 / dy dx 1 dx 2 f(\x 1 \, \x 2 \, |x 2 - « 

JO JO Jzi 

/■AT pN—qi 

6N dq! dq 2 f(\qx\,\q2\,\q2 + qi\), (36) 

Jo Jo 

provided that /(cci, x 2 , £3) is invariant under permutations of the arguments. The last 
equality is again justified in Appendix A In order to be able to sum up over fcj we use 
the following integral representations: 

1 



ioo+0 

T(h + k 2 + k 3 - 1/2) " 



i /-ioo+0 

— — / ds e s s -( fc i+ fc 2+fc3)+i/2 
2ni J_ ioo+0 



1 * 



+o 

rf/ ^2(fc 1+ fc 2 +fc 3 )-3_ ,3^ 

2(^ + ^ + ^-1) Jo 
The summation over k±, k 2 and k$ can be done easily now. All sums over ki have the 
form Ylk ) =o(~y) k ^r(l+i) ^ a = fot{v)i a = 0,1,2,3. The explicit expressions for f a are 
given by 



f Q (y)= h(y) 



y 



y(2 + y) y(4 + 2y + y 2 ) 

/2(2/)= -2(l + y)3/2 My) = ~ WTW 2 " ( } 

Using this notation and changing /3 by /3 = f3bt/ \/2 we can write P^ 3 ^ in a compact form 

wl dx L dyF A—-w 



g / d,/ W-fL JM, (39) 



where F 3 is defined as 



37r 3/2 6 2 



/iOO+U 
dse S y/s ^ Axiom x 
■ ioo +° ai, a ,3 



f ai (x 2 /s)f a2 (y 2 /s)f a3 CO* 6-1 + y e " 1 ) 2(e - 1) /s) • (40) 
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Note that Eg. (139 p is similar to the integral representation for P( 2 \ Eq. (jl8j) . 

6. Dynamical scaling exponent fi: the second order perturbation theory 
result 



To calculate the second order result for the dynamical scaling exponents /i one needs 
first to find 

n 3 (r) = lim (>> _ V 2 )) 2 ) . (41) 

For 

one can use the representation (1T91) allowing to integrate over (3 explicitly. To 
exploit a similar trick for P^ we scale the integration variables in Eq. (I39[) x — > (3~x, 
y f3~y. Then r-dependence get factorized for both terms in ( 141]) and the derivative 
with respect to In r can be now taken explicitly: 



<9IT 

T l-c 



J™dxj™dyF 3 (x^jT 1 ) - J 2 



(42) 



91nr 

with J defined in Eq.fl20D. 

The e-expansion of J is given by Eq. (I22|) ; its leading term is a constant of order 

of b. Therefore (— ) J 2 diverges in the limit e — > as 1/e. We show below that the 

first term in the brackets on the r.h.s. of Eq. p2|) also contains a divergent contribution 

of order 1/e , so that two divergent contributions cancel. This cancellation is another 

manifestation of the existence of the critical dynamical scaling. 

Now let us find an e-expansion for the double integral I = J °° dx f °° dyF 3 (x e_1 , y e ~ 1 )- 

i i 

To this end we change variables x by s 2 < 1 - e )x, y by s ^-^y, where x and y are com- 
plex. Then the integration over s in Eq. (14131) leads to the appearance of the inverse 
T-function (137|) . Next we deform the contour of the integration over x and y back to 
the real axis and get 

2 r (i^ - I) 7o x^V^^ + y) 1 - 

G(x, y, e) = ^ A aict2Ct3 g ai (x 2 ^) g a2 (y 2 ^) g a3 ((x + y) 2 ^ 1 )) , (43) 

"1,2,3 

where the functions g a are related to f a defined in Eq. (l3"g|) as g a (y) = f a (y)/y/y~- in 
order to find the e-expansion of the integral over x and y in the above equation, it is 
convenient to change the variables by q = x + y and z = x/(x + y). In terms of the new 
variables the integral, which we denote by Jo, takes the form: 

f 1/2 f°° da 1 
/o = 2 1 dz J -^ z^{l-z)^ G{qZA{1 - zU) - (44) 
To derive this equation we used the fact that G is a symmetric function, i.e. G(x, y) = 
G(y,x). The leading singular term of the e-expansion of Iq originates from z — > and 
can be extracted by integration by parts: 



z JO 1 



G(qz,q(l 



=1/2 
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1/2 d 
dzz e — 
dz 



U-e 



Q 



dq 

2-3e 



G(qz,q(l - z),e) 



(45) 



.4 



Thus the e-expansion of Iq has the form Jo = (2/e)L4 + e£> + 0(e 2 )]. The coefficient A 
is obtained by setting e = in all the terms in the brackets in Eq. (145 p and it equals to 

^G(0, g ,0) = |, (46) 
6 

where the integral is calculated explicitly using the definition of G (143j) . To calculate 
i? we collect the first order terms in e generated by all e-dependent contributions in 
Eq.fllH]) and find 

B = ---]n2 + R. (47) 

2 6 V 1 

The first two constants in the above formula are equal to integrals over q similar to one 
in Eq. pBl) . The constant R is given by the two-dimensional integral 



R 



1/2 l 

Z 



'l-z) 



dq 



G(qz,q(l-z),0) 



0.276, 



(48) 



7T 

" " ' 6 

which we were able to compute only numerically. The derived results for A and B along 
with the e-expansion of the inverse T-function in Eq. (1430 yields 



n 2 b 2 
~2e~ 



1 + e 3 + 7 + ln2 



6R 

71 



0(e 2 



Substituting this formula as well as the expression for J Eq. (I22j) into Eq. 



dU 3 2l2 f3R In 2 

— = 71 

dlllT \ 71 2 



-0.8196 2 . 



(49) 
we obtain 
(50) 



An alternative integral representation of this answer can be found in Eqs. (22-23) of 
Ref.[12]. We emphasize that the singular 1/e terms cancel giving a finite result in the 
limit e — > 0. Thus we have demonstrated the existence of the dynamical scaling with 
the accuracy of the sub-leading terms of the perturbation theory. 



7. Fractal dimension g^: the second order perturbation theory result 

To calculate the fractal dimension di one needs to deal with 



n 3 (A0=Hm (p (3) - ^(P (2) ) 2 ) • 



(51) 



Changing the integration variables x = /3 1 - e Nx, y = fi^-^Ny and (3 = N (3 in the 
integral representations fTl~8]) and (|39|) for and respectively and taking the 
derivative with respect to In N we find 



dn 3 

dlnN 



2eN 



2c 



3s-l 

d(3(3— 




2e-l 

d(3P^ 



dx 



dxFi 



dyFs 



x 1 e ' y 1 e 



x 



(52) 
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where we returned to the previous notation for x, y, and (3. Since we are interested in 
the limit e — > we can integrate over using the formula f l2~Tj) . In this way we obtain 
that q^ 3 n is very similar to the expression (H2|) for 
f)T\ f)T\ 

lim ^ = lim ^ + lim[2ei^(e)] - K 2 , (53) 

with 

ijTi(e) = / d/31n/3/3^ / dxF 3 — , , (54) 

Jo Jo \ x ((3 i-« - x) l ~ e J 

K 2 =2 [ dxF 2 (-) [ df3\nf3^-F 2 (f3). (55) 



jo \ X J Jo ft 2 

The r.h.s. of Eq. fl53|) does not contain divergent contributions of order 1/e, as 
divergences cancel out in ^ (see the previous Section) and lim e ^. [2e-R'i(e)] is finite 
(see Eq. (l59|) below). This fact demonstrates the existence of the spatial scaling 1/N d2 
with the accuracy of the sub-leading terms of the perturbation theory. The validity of 
the Chalker's ansatz implies that lim e ^o[2eA'i(e)] — K 2 = and this is what we show 
below. 

Let us first calculate K 2 . Using the explicit expression for F 2 ( |T8i) one can easily 
evaluate both integrals 

dxF 2 ( — | = — -^tL 
d/31n/?-F 2 (/3) = -/=(1 + 7 /2 + In 2), (56) 

hence we obtain 

J^ 2 = - 7r 2 6 2 (l+ 7 /2 + ln2). (57) 

For Ki we are interested only in the leading 1/e term of the e-expansion. In order to 
extract it , we first change the variable p by y = /3 — x in Eq. (1541) and find that 
K x {e) = -{l-efh with 

h = So ^ So dyHx + y)Fs {'x^^)- (58) 
The expression for I\ has a structure very similar to the structure of I defined below 
Eq.( )42|) and hence its e-expansion can be found in exactly the same way. Skipping the 
details of the calculation we present here only the final result 

lim[2eATi (e)] = -n 2 b 2 (l + In 2 + 7 /2) = K 2 . (59) 
t— >o 

Thus we conclude that the contributions of the two last terms in Eq. (I53[) cancel and 
lim^o = li m e->o^§7-- This equality not only proves the validity of the Chalker's 
ansatz in the second order perturbation theory, but it also provides the expressions for 
d 2 and \x: 

Tib 9 , 9 /3i? In2\ _,,q, , i 

d 2 = l-n= — - n 2 b 2 I — - —j + 0(b 3 ), (60) 

where R is defined in Eq. 
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8. Conclusion 

In the above calculations we have demonstrated by expansion in the parameter b <C 1 
that the power-law scaling Eq.flT]) and Eq.(j3]) holds true as soon as IniV ^> 1 and 
ln(Eo/u) ^> 1, Eq ~ b, even in the perturbative region where b IniV <C 1 and 
b \u{Eq/uj) <C 1. This statement is verified up to the second order in b <C 1- With 
the same accuracy we have shown that the exponents di and fi are connected by the 
scaling relation Eq.(jl]). Moreover we have found a term ~ (tt&) 2 in d<i (see Eq.( !60l) ) 
which appears to enter with an anomalously small coefficient 0.083 (71b) 2 . 

However in order to obtain all the above results we used the analogue of the 
dimensional regularization, replacing (b/(n — m)) 2 in Eq.Q by (b/\n — m]) 2 ^ 1 "^ and 
setting e — > at the end of calculation. This trick is well known in quantum field theory 
and it is working well for models whose renormalizability is proven. In other words, it 
works well, if it is known that the critical exponents (and the power-laws themselves) 
do not depend on the short-range details of the system (e.g. on the form of the function 
(\H nm \ 2 ) in Eq.Q which interpolates between the well defined limits at n = m and 
\n — m\ ^> 1). We would like to emphasize here that the renormalizability of the long- 
range model studied in this paper is not proven. That is why it may in principle happen 
that the results derived in the present work depend on the regularization scheme. We 
only know that this is not the case in the first order in b where all the integrals can 
be explicitly calculated using any other regularization. Whether or not the universality 
(independence on the interpolating function) holds in higher orders in b is an interesting 
open problem. 
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Appendix A. Integration over the "center of mass" coordinate 

Passing from discrete sums to integrals in Eq.( |T5l) and Eq.( l36l) . we replaced the 
integration over the "center of mass" coordinate y by multiplication by N. The aim of 
this Appendix is to justify that step. 

In calculating we deal with the following integral: 



1 r" rv 1 r n 

h(N) = -J dy J dxf(x) = -J dyF(y), 




(A.l) 



Our calculations show that the asymptotic behavior of F(y) is given by 
F(y) = c\ny + c + O{l/y), 



(A.2) 
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hence 

I 2 (N) = c(ln A - 1) + c + . . . = cln(JV/e) + c + . . . « F(N/e). (A.3) 

So that replacing (1/iV) dy by one in Eq. flA.ll) is equivalent in the asymptotic limit 
to replacing A by A/e. Now let us show that the same is true in calculation of P^(t). 
The relevant integral has now the following form: 

I 3 (N) = —J dy dq 1 I dq 2 f(\q 1 \,\q 2 \,\q 2 + q l \) = —j G(y). 

(A.4) 

According to our results the asymptotic behavior of G is given by 

G(y) = d 2 ln 2 y + d x \ny + d + 0(l/y), (A.5) 

then substituting this expansion into the definition of /3(A) we find 

I 3 (N) = d 2 [ln 2 A-21nA + 2]+d 1 [lnA-l] + d + ... = 

= d 2 ln 2 (A/e) + d x ln(A/e) + d + . . . « G(N/e). (A.6) 

Thus the integration over the "center of mass" can be taken into account in calculations 
of both p( 2 ) and P^ by scaling the system size. However since we are actually interested 
in calculation of the scaling exponents the scaling of the system size by a constant factor 
is irrelevant as it follows from Eq.(j9]). 

Appendix B. Scaling exponents in the orthogonal symmetry class 

The leading contribution to the virial expansion of the return probability in the 
orthogonal case can be obtained straightforwardly from the results of Ref . [16] : 

[k- n 

P£L(t) = - V £ ^^^Pnltllo (2M 2 ) . (B.l) 

Here Iq is the modified Bessel function. We have to calculate the double sum in Eq. (IB . 1 1) 

( 2 \ 

with logarithmic accuracy at bt 3> 1 and N 3> 1. Therefore the formula for P orth can 
be reduced to a single integral (cf. the unitary case): 

la.-^jf^g).,.-, ( B,) 

where / is a finite constant. The dominant contribution to the integral is governed by 
the region x < r where the asymptote of the Bessel function can be used: 

Uz > 1) 



Thus we can rewrite Eq. (IB.2l) with the logarithmic accuracy as follows 



min(7V,r) dx 



Inserting Eq. (IB.3l) into Eq.Q we find 

1-H = d 2 = V2b + 0(b 2 ), (B.4) 
which confirms the Chalker's ansatz up to the leading terms of the perturbation theory. 
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